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<^ ; Abstract 

\^ ' The very recent analysis by BaBar Collaboration indicates that the X(3872) may favor the 

^^ ■ quantum number J = 2 ^ rather than the previously assumed l"*"*". By pretending the 

?7c2(l-D) charmonium to be the X(3872), we study the parity-even radiative transition processes 
?7c2(l-D) — )• J/'0('0') + 7 within several phenomenological potential models. We take the ^Di ad- 
mixture in ?/^' into account, and consider the contributions from the magnetic dipole (Ml), electric 
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Qj ' quadrupole {E2), and magnetic octupole (M3) amplitudes. It is found that the ratio of the branch- 

f^, ' ing fractions of these two channels, as well as the absolute branching fraction of r/c2 — ?• '4''1j ^^^ 

in stark contradiction with the existing BaBar measurements. This may indicate that the 2 '" 
assignment for the X(3872) is highly problematic. 
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The rising of a dozen of new charmonium resonances in recent years, most of which are 
above the DD mass threshold, has greatly reinvigorated the field of hadron spectroscopy. 
Among them, the X(3872) particle, perhaps being the most intensively studied one from 
both theory and experiment, has occupied the central stage [if. The X(3872) was first dis- 
covered in 2003 by Belle Collaboration in the decay B'^ — )• J/ipn'^TT^K'^ |2[, subsequently 
confirmed by Babar J3|, as well as CDFII y^, DO |5| in inclusive pp collision experiments. 

Several unusual properties of this particle, i.e. its mass in extreme proximity to the D D 
threshold, the rather narrow width, and the large isospin violation seen in its decay pattern, 
have persuaded many authors to believe that, rather than being a conventional charmo- 
nium |6|, the X(3872) may be of exotic nature, e.g., a loosely-bound D'^-D molecule [Tf, 
or a diquark-antidiquark cluster |8|, or a hybrid |9|. 

Aside from its mass, the most important property about the X(3872) is its J^*" quantum 
number. The even C-parity of the X(3872) has been firmly established since the recent 
observations of its decay to J/ijjp^ [lO[ and to J/ip'y [ll[. By analyzing the decay angular 
distribution in the process X(3872) — )■ J/iIj7t^tt~, CDF Collaboration narrowed the pos- 
sibilities of the J^ down to 1''" or 2~ [12']. The former assignment seems more appealing 
from the theoretical perspective, which is particularly congenial to a S-wave D D molecule 
interpretation. Although the smoking gun from the experimental side is still not available 
yet, the l^"*" assignment for the X(3872) has already been tacitly accepted in most of the 
recent phenomenological works. 

However, very recently there comes a quite unexpected, and, perhaps disquieting, news 
from BaBar Collaboration. The latest analysis of the decay B — >■ J/ipuK by BaBar 
indicates that the P-wave orbital angular momentum for the J/ipu system is more favored 
than the S-wave, which implies that the X(3872) may favor J^"-^ = 2 '" instead of the 
universally-believed l^"*" [13|. Therefore, we are compelled to re-scrutinize the properties of 
the X(3872). If future experiments will confirm the result of [l3|, our perception on the 
nature of this particle would have to be profoundly changed. 

In light of the latest BaBar analysis [13|, the most natural candidate for X(3872) would 
be the ?7c2(l-D) meson. This D-wave spin-singlet (denoted by the spectroscopic symbol 
^D2) charmonium has been extensively studied in quark potential models for decades. Its 
predicted mass is scattered in the range 3760-3840 MeV [6|]. The width oir]c2 is believed to be 
narrow, since the decay into DD is forbidden by parity, and the energy conservation does not 
allow it to disintegrate into DD . It can only decay through the strong and electromagnetic 
transitions, exemplified by ric2 -^ Vc'^'^ ^^'^ Vc2 -^ hc'~f, as well as through the OZI-forbidden 
annihilation r]c2 -^ gg. Each of these processes is expected to have a partial width of a few 
hundred keV. 

One of the strongest objections to identifying the X(3872) with the 77^ 2 charmonium 



probably comes from the electromagnetic transitions ?7c2 — > J/i/j{ip') + 7 [1^. Such parity- 
conserving transitions fiip the quark spin and change the orbital angular momentum by 
two units, so there must be strong multipole suppression, and one would expect a rather 
small branching fraction for such decay processes. Thus, one is puzzled by the fact why 
the BaBar Collaborations were able to observe these radiative decay channels several years 
ago ll|, with only a limited statistics of the X(3872) samples? 



The motif of this paper is to quantify this objection, by presenting a detailed study for the 
radiative transition processes 77^2 — ^ J/ip{ijj') + 7. We will employ potential nonrelativistic 
QCD (pNRQCD) as a convenient calculational device, and work with several phenomeno- 
logical potential models, as well as take the ^Si — ^Di admixture effect for tp' into account. 



We will also identify the contributions from several different multipoles, e.g., the magnetic 
dipole (Ml), electric quadrupole (-E'2), and magnetic octupole (M3) amplitudes. 

Our key finding is that, no matter the mixing effect for ip' is taken into account or not, 
the predicted branching ratio for 77^2 -^ ip' + 'y is orders of magnitude smaller than the lower 



bound that can be inferred from the BaBar measurement [15l]. Moreover, the ratio of the 
branching fraction of 77^2 -^ ip' + '^ to that of ?7c2 — >■ J/ip + 7 is orders of magnitude smaller 
than the corresponding BaBar measurement for the X particle [l5|. These qualitative 
conclusions do not vary with specific potential model. We thus tend to conclude that, if the 
BaBar experiment is correct, the 7]c2 assignment for the X(3872) particle would become 
highly unlikely. We hope that future higher-statistics experiments will help to clarify the 
situation. 

Before launching into the calculation, we first recall some related experimental facts. 
BaBar has recently measured the following products of two branching fractions [l5|: 

B[B^ -^ X{3872)K^] B[X{3872) -^ J/tP + 7] = (2.8 ± 0.8 ± 0.1) x 10-^ (la) 
B[B^ -^ X{3872)K^] B[X{3872) ^ ^' + 7] = (9.5 ± 2.7 ± 0.6) x 10-^ (lb) 

with 3.60" and 3.5a significance, respectively. Therefore one can deduce the ratio of the 
branching fraction of X — )■ ip''-/ to that of X — )■ J/ip'j: 



i3[X(3872) -> J/i/j + 7] 

This measurement is in serious conflict with the predictions made from some specific D D - 
molecule models [7|, but the calculational framework underlying those models may be ques- 
tionable. Interestingly, this ratio seems roughly compatible with the canonical Xci(2-P) 
interpretation of the X(3872), which decays to J/ip{ip') + 7 through the dominant electric 
dipole (El) transition |6l. Il6l-ll8| ^. 



Very recently BELLE Collaborations have also analyzed these two radiative transition 
channels [l9[. Their preliminary results are 

B[B^ -^ X{3872)K^] -B[X(3872) ^ J/tP + 7] = {1.78toil ± 0.12) x 10"^ (3a) 
B[B^ -^ X{3872)K^] -B[X(3872) ^ V'' + 7] < 3.4 x 10"^ (3b) 



Their first measurement is consistent with BaBar's result, flTa|) . But Belle has not 
observed any i?^ — )■ ip''yK'^ signals. Consequently, Belle is only able to place an upper 
bound on the ratio of these two branching fractions: 

S[X(3872)^^' + 7] ^2.1, (4) 



i3[X(3872)^ J/^ + 7] 



at 90% confidence level. 

Thus far, our knowledge is only limited to the product of two branching fractions, and it 
will be certainly useful to know the absolute branching fraction of X — ;■ J/ilj{ip')'y. Using a 



^ It seems that the ^Si—^Di mixing effect has not been incorporated in most phenomenological analysis of 

Xci(2P)^V'7- 



missing mass technique, BaBar was able to set an upper bound for the absolute branching 
fraction of B^ -^ X(3872)Js:± some time ago |20|: 

B[B^ -^ X{3872)K^] < 3.2 x 10-^ (5) 

at 90% confidence level. 

Combining the results of ([1]), ([3]) together with ([5]), we are able to place some lower 
bounds on the absolute branching fractions of the X(3872) decays to J/ip{ip') + 7: 

B[X{3872) -^ J/i) + 7] > 5.9 x 10^^ i3[X(3872) ^ ^' + 7] > 1.9 x 10"^ Babar 

(6a) 
S[X(3872) -^ J/i) + 7] > 3.8 X 10"^ Belle (6b) 

We are unable to establish a meaningful inequality for X — )■ ^' + 7 from the latest Belle 
data. 

Our central task then is to explicitly examine whether the decays ?7c2 — )■ J/tpiip') + 7 are 
consistent with the experimental constraints listed in ([2]), (^, and (|6]). 

In passing, we note that there have already existed some theoretical studies on the decay 
Vc2 -^ J/ ^{4'') + 7 by Sebastian and his coworkers 2l|-|23|. These authors worked in the 



traditional first-quantized quantum-mechanical framework. Notably, some of serious incon- 
sistencies seem to exist amongst a sequence of their papers ^. Therefore we feel it timely 
and obligatory to conduct an independent investigation on these transition processes. 

A modern effective-field-theory framework of dealing with the single-photon transition 



in quarkonium has recently been put forward in Ref. [2J], by coupling the potential non- 
relativistic QCD (pNRQCD) |25| with the electromagnetism. In [2J], the Ml transition 
process with A/ = 0, exemplified by J/ip — )■ rj^'y, has been systematically analyzed. In this 
work, we will utilize the same formalism to tackle the A/ = 2 magnetic transition process 
Vc2 -^ J / i^ii^')! ■ Our final results turn out to significantly differ from those in |2lM23l . 



To account for the radiative transitions, it is convenient to promote the gauge group of 
NRQCD to S'f/(3)c X f/(l)cm- One then matches this NRQCD action onto an even lower 
energy effective field theory, pNRQCD, by further integrating out the quantum fluctuations 
of virtuality of order im?v'^. pNRQCD is an ideal formalism to tackle radiative transitions, 
because the active degrees of freedom, i.e., the dynamical gluons and the emitted photon, 
are both ultrasoft [k^ ~ mv"^) and can be treated on an equal footing. Since they possess a 
wavelength much longer than the typical quarkonium radius, the corresponding fields need 
be multipole-expanded. By this way one can elegantly implement the multipole expansion, 
the standard treatment of the electromagnetic transitions in quantum mechanics textbooks. 

Out starting point is the following pNRQCD Lagrangian density |24l |: 



/^pNRQCD = [ d\ Tr h^ (ido + ^ + ^- V^''\r)] S I + Aight + C.^^nacD , (7) 

where 5* represents the field for a composite system made of the heavy quark and the heavy 
antiquark, depending on the center-of-mass coordinate, R, and the relative coordinate, r. 5* 
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For example, the partial width of i]c2 — !• J/i^ + 7 was originally predicted to be 2.13 keV 2ll| . This 



prediction later shifted to an abnormally large value, 62.6 kcV [23. In their last publication [23[, this 
prediction has shrunk to 0.699 keV. 



+ ^T£^{S^-■(r■V.fB<=-}S + £^rrr {St,.-rx(fxBn}S 



is a 2 X 2 matrix in spinor space and a singlet under color and t/(l)cm transformations. The 
Vg denotes the static singlet QQ potential. The first term in ([7]), describing the dynamics 
of the QQ pair dictated by the strong interaction, resembles very much the traditional 
potential model. £iight represents the lagrangian for light degrees of freedom, e.g., light 
quarks, gluons as well as photons. The particularly relevant to this work is the third term, 
£^pNRQCD5 which depicts the spin-dependent interaction between the QQ pair and a photon: 

^C^PNRQCD = jd'rTil ee^Str ■ E^-S + ^^^ {S^, o" ■ B'^-} S 

em 

'f^ {St, ^ . (r ■ VnfB-} S + |^ -^^^'^' 
lorn ^ ^ 8m^ 

cm cm 

-T6^ [St,..HVHX,E^-]]S-^ [St,..HV.x,(r.V^)E<=-]]S 

cm /-i^^^zDzD 

where the trace is over the spin indices. The W^ and B*^™ signify the electric and magnetic 
field strengths, which depend only on R. The ubiquitous occurrence of the Pauli matrix 
cr (the first operator, responsible for the El transition, is an exception), signals that the 
concerned radiative transition is of the spin- dependent type, i.e., with quark spin flipped. 
The appearance of r is a consequence of multipole-expanding the electricmagnetic field, 
while the f = r/r represents the unit radial vector. Notice that ([S]) is manifestly gauge 
invariant. 

The coefficients c^, c|™, c^^g? ^^"^ ^"^ in ([8]) are the matching coefficients that were 
directly inherited from the NRQCD with enlarged gauge group. They satisfy some exact 
relations dictated by reparametrization (or Poincare) invariance (26|: 

em _ 9„cm i em _ em em _ i em _ cm i /q\ 

Cg — Z,Lp i, (^W12 — ^Wl ^W2 — -■-! S'P ~ ^F -■-• \^J 



All these matching coefficients are known at the one loop level [26|. In particular, we have 

4 as 



c 



1 + «^q" = l + o77^ + C?(«s), (10) 
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where Kq^ is usually dubbed the anomalous magnetic moment of the heavy quark. The one- 
loop perturbative contribution to Kq^ in ( ITOl) seems insignificant, less than 10% for charm. 
However, it is often conjectured that the magnetic moment of the bound quark may receive a 
large contribution due to some nonperturbative mechanism. For simplicity we will suppress 
such a possibility (see |2J| for a discussion on whether such a contribution can naturally 
emerge from the first principle of QCD). Rather we will content ourselves with taking kq as 
its short-distance value. 



The cr ■ r X (f X B''™) term in ([S]) is also worth some remarks. As discussed in [24 
unlike the remaining operators listed in (jH]), this spin-dependent transition operator receives 
some nontrivial corrections when descending from NRQCD to pNRQCD. This operator is 
intimately related to the quark spin-orbit potential, whose coefficient is also protected by 
the reparametrization (or Poincare) invariance (Gromes relation) |27j. 

Prior to giving concrete expression for 77^2 — )■ J/tp{'ip')'y, we note that there have been 
some phenomenological inclinations that ijj' may not be a pure 5'-wave vector charmonium. 



Rather it may have some admixture of Di component [28|: 



li^') = cos(p\2^Si)-sm(P\VDi), (11) 

usually with taken to be around 12°. 

To compute the radiative transition process in pNRQCD, it is necessary to know the 
corresponding vacuum-to-quarkonium pNRQCD matrix elements. To the purpose of this 
work, we need know 

{0\S{R,r)\n'S^{P,\)) = ^ i?no(r) ^-^^f^^ e'^'^ (12a) 



(0|5(R,r)|n^Pi(P,A)) = /J^ni(r) ^^^^^i|^ e^^'^, (12b) 

{0\SiR,v)\n'D2{P,X)) = J^Rn2ir) "^^;; e'^'^, (12c) 

V ovr V2 

(0|5(R,r)|n^A(P,A)) = -^i?..(r) ^^^^-^^-^~^^^^^-^ e-^-^. (12d) 



The first two entries have been constructed in [2J], and the last two are new. Here Rni{r) 
denotes the radial Schrodinger wave function of a quarkonium. The symbol e(A) repre- 
sents the polarization vector of any J = 1 quarkonium state, satisfying the orthogonality 
condition e*(A) ■ e(A') = Sw, whereas the symbol h^^{X) denotes the polarization tensor 
of the n^D2 state, which is symmetric and traceless, obeying the orthogonality condition 
tr [/i*(A)/i(A')] = 6\y. The overall normalization factors for each entity are chosen such that 
all the quarkonium states are nonrelativistically normalized. 

With the knowledge of ([8]) and (fT2l) . it is then a straightforward exercise to calculate the 
transition amplitude for ^D2 — > ^5*1 + 7 and ^D2 — ?■ ^Di + 7. The latter is a regular Ml 
transition between the D-wave spin singlet and triplet, which is included here because it can 
contribute to ?7c2 -^ ip' + 'j through the mixing mechanism. After completing the angular 
integration, the desired results are 



ff 2 r e* ■ k X £* k^W^k^ 
M[^D2^'S^ + ^] = ^^Y^l ^2 cTJi + k'h^' {^* ^ e^y [cT -l)J2 

+ e*'h^^ (k X ey {J2 + Ja- c^Js) }, (13a) 

M[^D2 -^ ^D, + 7] = ^7p ^"h'' (k X e;Y c^Jo, (13b) 

where k and e* represent the three-momentum and polarization vector of the emitted 
photon, respectively, and e* stands for that of the outgoing J = 1 charmonium. The 



involved dimensionless overlap integrals are given by 



oo 



Jo = / drR3D,{r)RiDAr)r\ (14a) 

Ji = —i- / drR^ssAr)Rw,{r)r\ (14b) 
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J2 = ^ drKss,ir)RwAr)r\ (14c) 

^"^Q Jo 

1 f°° / R' (r)\ 

■J3 = --T dr (r:,sM) -^^)R^D.ir)r^ (14d) 

^Q Jo \ r / 

J4 = 7^ / rfri?„35,(r)yf (r)i?i,5,(r)r3. (14e) 



To facilitate the comparison with the expressions in Ref. |21|-|23|, we have introduced the 
same Jj (z = 0, . . . ,4) as theirs, except a different normalization factor for Ji is adopted. 

Some remarks on deriving the amplitude (I13ap are in order. Unlike in the case of an 
ordinary Ml transition process with A/ = 0, not all of the operators in (|H]) can make a 
nonvanishing contribution for such a A/ = 2 transition process. For example, the leading 
magnetic- dip ole operator, the cr ■ B°™ term, and the cr ■ [— iVijX, E''™] term, as well as the 
(T ■ B°™V^ term, simply cannot connect a D-wave state to the 5'-wave state owing to their 
spherically-symmetrical nature. It turns out that only four terms in (|8]) can directly make 
a nonzero contribution. 

Equation (I13ap has also embedded an interesting piece of contribution, the so-called final- 
state recoil correction |29|. As a manifestation of the relativistic effect to a moving ^Si state, 
its wave function can develop a nonvanishing overlap with a spin-singlet P-wave component, 
therefore the transition ^D2 — t- ^Si can be effective realized through a El transition from 
the parent to this small ^Pi component. Some subtle gauge-invariance issue related to this 
Lorentz boost effect has been clarified (2J]. Although the recoil correction and the spin- 
orbit-potential-related term are not separately invariant under the f/(l)em transformation, 
their sum is. We have explicitly verified that, for the reaction ^D2 — )■ ^Si + 7, the sum of 
the recoil correction and the correction induced by the spin-orbit-potential-related operator, 
is indeed independent of the redefinition of the pNRQCD field S by implementing a f/(l)em 
gauge link. 

For the decay ^D2 — ?■ ^Di + 7, which is an ordinary allowed Ml transition, we only 
consider the leading contribution and not include the relativistic correction in (I13b|) . 



Squaring the amplitudes in ( !T3|) . and summing over all possible polarizations, it is easy 
to obtain the spin-averaged partial width r[?7c2 — ;■ J/ipi'ip') + 7]- Nevertheless, (fT3|l encodes 
much richer polarization information. In this work, we would like to proceed to extract the 
helicity amplitude |30|, |3l| and the multipole amplitude |32|, |33| from this equation. These 



two types of amplitudes can in principle be extracted experimentally. It is worth noting that, 
CLEO-c experiment has recently extracted the higher-order multipole amplitudes associated 
the cascade decay process ip' — ?■ Xci,27 -^ J/"^ + 77, by performing a maximum likelihood 
fit of the joint angular distributions of the two photons |34| . 

For the decay ?7c2(i^) — > "^(z^) + 7(-^)) we signify the projection of the angular momentum 
of ?7c2 along the moving direction of the 7 by u, and denote the helicities of ip and 7 by /i and 
A. With this specific choice of the quantization axis, we have v = \ — ^. There are in total 
three independent helicity amplitudes A^^\, and the other three can be related by parity 



invariance. We introduce the shorthand A^, for A^ a with A fixed to be +1. Substituting the 
exphcit representation of polarization tensors in (TT3|) . it is easy to obtain: 



^0 = A, 



Ai 
A2 



Ao,i = 



-A-i = 

-A. I 



2cf^Ji - (24^- - 1) J2 + c^'^pJs -Ja + SV2cf^Jo sin 
-c7 J2 + cX J3 -Ja + 3V2cTJo sin (j) 



pp 



V2 
J2 + CIZJ3 -Ja + 3^24"^ Jo si 



-p'p 



smi 



(15a) 

(15b) 
(15c) 



where is the ^Si — ^Di mixing angle, which equals for J/ip^ and may be put 12° for 
ijj' on phenomenological ground. In deriving the helicity amplitudes for ip' + 7, we have 
approximated cos0 = 0.978 ~ 1 for simplicity. 

From (TT5|) . it is straightforward to obtain the spin-averaged partial width: 



T[7^e2^j/^m+i] 






A |2 _L I 4 |2 I I /I |2 
/IqI + |/1i| + |/l2| 



(16) 



According to [32|, |33|, the helicity amplitudes Ay {u = 0,1,2) are connected to the 



multipole amplitudes aj through the following orthogonal transformation ^■. 



A. = Y. 



2J^ 



2 J,.. + 1 



aj-y {JfA'A,'^- MJvc2^'^), 



(17) 



where the Condon-Shortley notation for the Clebsch-Gordan coefficients, 
{ji,7ni;J2,m2\J,M), has been adopted [35|. J^^^ = 2 is the spin of the rjc2 meson. 
J^ represents the angular momentum carried by the photon, obeying 1 < J^ < Jr,^2 + 1; ^^^ 
ai, 02 and a^ correspond to the Ml, E2, and M3 multipole amplitudes, respectively. 
The inverse transformation of flT7|) can be readily obtained: 
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From ( ITSl) . we can readily deduce these three multipole amplitudes: 

2cf'Ji - 5(1 + cf')J2 + lOc^^Ja - 10 J4 + 30^24°^ Jo sin, 



ai 
02 
03 



2V15 



2c5?^Ji-3(c|"^-l)J2 



2^3 



2crJi 



(18) 



(19a) 
(19b) 
(19c) 



Note that this transformation matrix is identical to the corresponding one in the Xc2 — > J /^ + 7 pro- 
cess 34|. But for such a i?l-dominant process, ai, 02, 03 should be identified with the electric dipole (-E1), 



magnetic quadrupole (A/2), and electric octupole {EZ) amplitudes, respectively. 



As expected, the S — D-wave mixing effect can only enter into the Ml contribution. Note 
the M3 amphtude solely receives the contribution from the cr ■ (r ■ Vj^;)^B°™ operator in 
dH]), which comes from multipole-expanding the leading magnetic transition operator to the 
second order in r. 

We are now in a position to compare our expressions for the three multipole amplitudes 



f ll9p with those reported in |22l . |23| . Somewhat surprisingly, it seems that our expressions 



disagree with theirs for each multipole, except for the Jq and Ji pieces. The authors of 



22l . |23| have not included the anomalous magnetic moment of the c quark, and in their 
formulas, all the overlap integrals, Ji through J4, contribute to the Ml and E2 amplitudes. 
Taking kq = in f lT9|) . we find that J3 would disappear from the Ml amplitude, and the 
E2 amplitude will only depend on Ji, which are in diametric contradiction to the equations 
in 



22, 23 



We can express the spin- averaged transition width for 77^2 -^ J /ip{ip') + 7 in terms of 
three multipole amplitudes: 

2ae^ |k|3 / \ 

r[r/,2^JM^')+7] = 75^2 [\<^i? + \<^2? + \a,?y (20) 

which is equivalent to flT6|l . since the transformation flT7|) preserves the norm. 

To make concrete predictions for the transitions ?7c2 — ?■ J/ipiip') + 7, we need know the 
radial wave functions of each involved charmonium from some phenomenological potential 
models. To ensure that our conclusion not to rest heavily upon one specific model, we 
choose to study with five different potential models. All of them only differ in the way of 
parameterizing the static potential V^ . Sp ecifically, the potentials we choose are Cornell 
type 36|, the Buchmiiller-Tye (BT) type_|38|, NR potential by Barnes et al. [37|, the screened 
confinement potential by Li and Chao [l7|, and the potential proposed by Fulcher 39|. We 
solve the Schrodinger equation with these potentials numerically, with the input parameters 
taken from the aforementioned papers. 

In Table [H we tabulate the predictions from the various models for the overlap integrals 
Jj (i = 1, ■ ■ ■ ,4). When evaluating these integrals defined in (1141) . we have used mc = 1.5 
GeV, and determined the photon momentum by physical kinematics, i.e., we assume ?7c2 
with a mass of 3872 MeV and use the physical masses of J /ip and ip' as input. Thus we 
obtain |k| = 698 MeV ^ for ?7c2 -^ J /'4^1 and 181 MeV for 77^2 — )■ ip'"^. From Table HI it is 
reassuring that these overlap integrals are not very sensitive to the different models. On 
the other hand, the overlap integral Jo can be safely put to be unity, since to the intended 
accuracy, both the ?7c2(l-D) and the ^(1-D) have degenerate radial wave function. 

In Table UTt we give the explicit predictions for the decay rjc2 -^ J/ip{ip') + 7, such as the 
partial width, the helicity amplitudes and the (normalized) multipole amplitudes for each 
decay channel within various potential models. We have taken a = 1/137 and the electric 
charge of the charm quark, Cc = 2/3. We have used the one-loop perturbative value for the 
anomalous magnetic momentum of the c, as indicated in fITU]) . If as{mc) = 0.35 is used, we 
then get Hc = 0.074. Detailed numerical checks reveal that all the predictions in Table |TT] 
only change modestly if we set Kc = 0. 



^ One may argue that, the 700 MeV photon in r]c2 -^ J/i^l is a httle bit too energetic, and it is perhaps 
more appropriate to count |k| '^ mv in this case. This may cast some doubt on the vahdity of multi- 
pole expansion, the underlying basis of this work. It might be illuminating to apply the hard- scattering 



mechanism, recently developed in 40|, to reinvestigate this channel 



From Table [TTl one finds that the predicted partial width for r]c2 — ;• J/'ip'y from various 
potential models ranges between 3.11 and 4.78 keV, that for t]c2 -^ 'ip''~f without considering 
^Si—^Di mixing ranges from 0.017 to 0.029 keV, and that for ric2 -^ ip''~f incorporating the 
mixing effect (0 = 12°) ranges from 0.49 to 0.56 keV. This seems to indicate that different 
potential models give rise to reasonably consistent predictions for each channel. 

TABLE I: The overlap integrals Jj for the electromagnetic transition ?/c2 — ^ J/'^i'^')^ in various 
potential models. Since Jq = 1, we have not listed its value. 



^^^ Potential 
^^^lodels 

Ji ^^ 


Cornell [36] 


Screened [17] 


NR [37] 


BT [38] 


Fulcher [39] 




J/i; ij' 


J/i; i:' 


J/i) V' 


J/i: ^' 


J/i: ij' 


Ji 


-0.600 0.123 


-0.710 0.180 


-0.728 0.161 


-0.757 0.153 


-0.763 0.147 


J2 


-0.376 0.051 


-0.347 0.063 


-0.365 0.056 


-0.383 0.048 


-0.390 0.044 


Js 


-0.304 -0.256 


-0.227 -0.160 


-0.242 -0.191 


-0.245 -0.212 


-0.249 -0.225 


J4 


0.136 -0.243 


0.121 -0.218 


0.128 -0.231 


0.144 -0.244 


0.173 -0.291 



TABLE II: The predictions of rjc2 — )• J/i^O;') + 7 from various potential models. The mixing 
angle i> has been taken for both 0° and 12° for ^'. We have taken a = 1/137, and Kc = 0.074 by 
using as{mc) = 0.35. In addition to the spin-averaged partial width, the helicity amplitudes and 
the (normalized) multipole amplitudes for each decay channel have also been given. 



Potential Models 




<P 


^0 Ai A2 ai 


02 


03 


a2/ai\ 


ks/ail 


Width (keV) 




J/i; 


~ 


-0.39 0.19 0.22 0.31 


-0.66 


-0.68 


2.15 


2.21 


3.11 


Cornell 


i;' 


0° 


0.17 0.12 0.17 0.93 


0.26 


0.25 


0.28 


0.27 


0.017 




12° 


0.56 0.79 1.12 0.99 


0.05 


0.05 


0.05 


0.05 


0.50 




J/i; 


- 


-0.50 0.18 0.21 0.19 


-0.70 


-0.70 


3.72 


3.70 


4.22 


Screened 


i;' 


0° 


0.21 0.09 0.14 0.85 


0.38 


0.37 


0.45 


0.44 


0.017 




12° 


0.60 0.76 1.09 0.99 


0.07 


0.07 


0.07 


0.07 


0.49 




J/i: 


- 


-0.50 0.19 0.22 0.20 


-0.69 


-0.69 


3.49 


3.49 


4.45 


NR 


i;' 


0° 


0.20 0.11 0.16 0.89 


0.33 


0.32 


0.38 


0.36 


0.018 




12° 


0.59 0.78 1.11 0.99 


0.06 


0.06 


0.06 


0.06 


0.50 




J/i: 


- 


-0.53 0.20 0.22 0.18 


-0.70 


-0.69 


3.76 


3.76 


4.78 


BT 


i;' 


0° 


0.20 0.12 0.18 0.91 


0.30 


0.29 


0.33 


0.31 


0.020 




12° 


0.59 0.79 1.13 0.99 


0.06 


0.06 


0.06 


0.06 


0.51 




J/i: 


- 


-0.54 0.18 0.20 0.14 


-0.70 


-0.70 


5.16 


5.16 


4.77 


Fulcher 


i;' 


0° 


0.22 0.16 0.23 0.94 


0.24 


0.23 


0.26 


0.24 


0.029 




12° 


0.60 0.83 1.18 0.99 


0.05 


0.05 


0.05 


0.05 


0.56 



In contrast to the minor role played by retaining Kc, the effect of ^^i — ^Di mixing is 
enormous for the decay t]c2 -^ V^'7, as can be clearly seen from Table [TTl The predicted partial 
widths with = 12° are about 25 times larger than those without including the mixing for 
i:'. This can be attributed to the fact that the transition ^D2 — )■ '^-Di7 is an allowed Ml 
transition, with the overlap integral Jo = 1. However for the transition ^D2 — ?■ 2'^ 5*1 -|- 7, 
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due to the existence of a node in the 25* radial wave function, the corresponding overlap 
integrals are generally small. Therefore, even with a relatively small angle, the S—D-wave 
mixing can already play a very important role. 

In Table mi we also list the magnitudes of the Ml, E2, M3 amplitudes for each transition 
process. For clarity, we have renormalized a, (i = 1, 2, 3) in (TTUj) such that the new multipole 
amplitudes satisfy |aip + |a2p + jasp = 1. It is a very interesting observation that E2 and 
M3 amplitudes are almost identical in all the decay channels. For 77^2 — !■ J/'ip'y, these two 
higher-order multipoles are even more important in magnitude than the Ml amplitude! This 
pattern no longer holds true for the ?7c2 — ?■ ^'7 channel. Nevertheless, it is also interesting 
to note that, for r]c2 — )■ V^'7, the ratios |a2/ai| and jas/ail decrease significantly after the 
mixing effect is included. This can also be understood by the fact that, since the allowed Ml 
transition ^D2 — )■ '^-Di7 makes more pronounced contribution than ^D2 — )■ ^«S'i7, including 
the mixing effect thus can greatly amplify the importance of the Ml amplitude relative to 
other two multipoles. 

Surveying Table [ITl one may be able to place the following saturation bound for the ratio 
of the two branching fractions: 

g[,.2^V;- + 7] < ^ ,2° (21a) 

B[ric2 ^ J/V + 7] " 

^f-^,t;+"U 6.1x10-3, = 0° (21b) 

where the maximum ratio in first equation comes from the prediction in Cornell potential, 
and that in second equation is obtained from Fulcher's potential model. These results seem 
to seriously conflict with the corresponding BaBar measurements, ([2]), no matter the mixing 
effect for ip' is taken into account or not. This seems to be a strong evidence to disfavor the 
?7c2 assignment for the X(3872). At this stage, these upper limits seem still compatible, but 
orders of magnitude less, with the preliminary Belle upper bound on this ratio, (jl]). 

It is certainly interesting to examine whether the absolute branching fractions of r]c2 -^ 
J/ip{ip')'-f are compatible with the B factory measurements or not. To realize this goal, one 
first needs know the full width of the X(3872). Belle has set an upper limit for the total 
width of X(3872): F^ < 2.3 MeV at 90% confidence level Q. The 2008 PDG compilation 



has estimated the total width of the X(3872) to be F^ = S.Otl^ MeV [35 . 

We hope to establish an upper bound for the predicted branching ratios of ric2 -^ 
J/ip{ip')l- To this purpose, we have taken a conservative attitude, by assuming F^ = 1.3 
MeV, the lower end of the PDG estimate. Such a value by itself is consistent with the full 
width of the ?7c2. From Table UTt one can see that, among all the potential models analyzed 
in this work, the BT potential model gives the largest prediction to the partial width of 
Vc2 — ^ J/i'l, 4.78 keV. We thus estimate 

B[ric2 -> J/i^ + 7] < 3.7 X 10~3_ (22) 

It is interesting to contrast this theoretical upper bound with those inequalities inferred 
from the two B factory experiments, (E]). Equation (122 p seems to mildly conflict with 
the experimental lower bound from either BaBar or BELLE. Due to large theoretical 
uncertainties in F[?7c2 — ?■ J/ip + 7] and Tx, we are inclined not to make strong claim just 
based on the analysis for this channel. 

The decay ?7c2 — )■ '^'7 poses a much more stringent challenge to the data. Among all the 
potential models, Fulcher's potential model yields the largest prediction to the partial width 
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of ?7c2 -^ "^'7, with or without including the mixing effect. Taking the corresponding partial 
widths from Table [TTl we may place the following inequalities: 



B[vc2 -^^P' + l]< 4.3 X 10-^ 
B['nc2 ^f + -f] <2.2x 10"^ 



:12° 

0° 



(23a) 
(23b) 



Comparing these theoretical upper bounds with the lower bound imposed by the BaBar 
measurement, 1.9 x 10^^, which is given in ([6]), one observes the glaring discrepancies, no 
matter the mixing effect for ip' is included or not. If BaBar result can be trusted, this may 
also be viewed as a strong evidence against assigning the X(3872) as the ?7c2. 



TABLE III: The predictions of the overlap integral (r) and the partial width for the electric 
transition r]c2 — )• /ic + 7 in various potential models. The mass of he is taken to be 3525 MeV j35l ]. 



^^^^ Potential 
^^^odels 


Cornell 


Screened 


NR 


BT 


Pulcher 


k 

(MeV) 


(r)2i r 

(GeV"^) (keV) 


(r)2i r 

(GeV^^) (keV) 


{r)2i r 

(GeV-^) (keV) 


(r)2i r 

(GeV~^) (keV) 


{r)2i r 

(GeV~^) (keV) 


331 


3.06 587 


3.54 786 


3.39 720 


3.37 712 


3.32 692 



For the sake of completeness, finally we would like to contrast the ric2 ~^ J/4^{'^') + 7 
processes with the dominant electromagnetic transition of ?7c2, the parity-changing electric 
transition ?7c2 — ^ /^c + 7- With the knowledge of the leading electric-dipole operator in ([8]) 
and the vaccum-to-/ic matrix element in fll2bp . one can readily reproduce the well-known 
El transition formula: 






15 



(r) 



21 



/ dr Rip^{r)RiD^{r)r^. 
Jo 



(24a) 
(24b) 



The contributions from the higher-order multipoles, M2 and £'3, as well as the relativistic 
corrections, neither of which are expected to be significant, have been neglected for simplicity. 
From Table UTTt one can see that the predicted partial width for ric2 — ?■ /ic + 7 in five potential 
models ranges from 600 to 800 keV ^. As expected, this transition rate is several orders of 
magnitude greater than that of ?7c2 — > J/ip{ip') + 7. 

In summary, we have presented a detailed analysis to the radiative transitions r]c2 -^ 
J/ip{ip') + 7 within various phenomenological potential models, employing the pNRQCD 
formalism as an elegant and efficient calculational framework. The major motivation of 
this study is to examine whether identifying the X(3872) with the ?7c2(l-D) charmonium is 
compatible with the B factory measurements of the radiative decay X(3872) — >■ J /ip{%p') -I-7. 
Our study reveals such an assignment would cause severe contradictions with the available 
BaBar measurements, either for the absolute branching fraction of X — ?■ ^' + 7, or the 
ratio of this branching fraction to that of X — )■ J /ip + 7. As a result, if the BaBar 



^ The El transition rates tabulated in Table IIIII seem to be somewhat larger than the often quoted 464 
keV 
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measurements [15| are trustable, the possibility for the X(3872) to be identified with the rjc2 
then becomes rather low. In our opinion, the X('3872) is most likely to carry J^*^ = l"'"^. 

A nuisance is that the preliminary Belle results [l9| on the decay X(3872) — !■ ip'j seem 
to confiict with the BaBar measurement 15|. It is important and urgent to resolve this 



experimental discrepancy, in order to unambiguously unravel the nature of the X(3872). 
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